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ABSTRACT 

We present mean reddenings toward 3525 RRO Lyrae stars from the Galactic 
bulge fields of the MACHO Survey. These reddenings are determined using 
the color at minimum V^-band light of the RRO Lyrae stars themselves and are 
found to be in general agreement with extinction estimates at the same location 
obtained from other methods. Using 3256 stars located in the Galactic Bulge, we 
derive the selective extinction coefficient Rv,vr = Ay / E(V — R) = 4.3±0.2. This 
value is what is expected for a standard extinction law with Rv,bv — 3.1 ± 0.3. 



Subject headings: dust, extinction — Galaxy: center — stars: statistics — sur- 
veys 
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1. Introduction 

Studies of stars in the central Galactic region are essential for an understanding of the 
Milky Way. For example, the questions of the structure and extent of the Galactic bar, the 
star formation history, and the distance to the center of our Galaxy can all be answered by 
analyzing various stellar tracers. Knowledge of the Galactic extinction is mandatory to begin 
investigations of these stars in a quantitative fashion. A good understanding of extinction is 
particularly important for the direction toward the Galactic Bulge and within the Galactic 
plane, where extinction is very significant. 

Extinction Ax changes as a function of wavelength A. The shape of the e xtinction curve 



i s cha racterized by the quantity Ry, the ratio of total to selective extinction; ICardelli et al. 



( 119891 ) have shown that the family of curves characterized by this single parameter provides 
a good fitting approximation to most extinction curves. Thus, if the value of Ry can be 
determined (e.g., from optical and IR photometry), then one can approximately find the 
properties of the entire UV-IR extinction curve. 

There are a few suggestions that the coeffi cient of selective extinction in the direction of 



the Galactic Bulge may be non-standard (e.g. iPopowskil 120001 ) . Anomalous extinction, i.e., 
extinction described by smaller Rv,bv or Rv,vr ratio than the standard value of Rv,bv — 
3.1, can significantly affect analysis of the stellar populations in the Galactic Bulge. Here, 
Rv,bv is the coefficient of selective extinction AyjEiB — V), with reference to the B and 
V filters, and Rv,vr = Ay/E(V — R), alluding to the V and R filters. Especially with the 
onset of a new generation of large telescopes, this problem becomes increasingly pressing and 
important, and one that can be addressed with this study. Here, optical V and i?-band data 
for 3256 RRO Lyrae stard3 from the MACHO microlensing experiment are used to find both 
the coefficient of selective extinction and to compute the absolute reddening near the Galactic 
bulge/bar. The value of Ry t BV is known to vary in certain directions, and reddening laws 
are normally derived from differential measurements of reddened and unreddened luminous 
early-type stars. However, as such stars are rare in the Galactic Bulge, the normal derivation 
of Ry t BV can n °t be carried out. The Ry t sv values determined from the RR Lyrae stars are 
compared with other reddening relations. 

The data used in this analysis are based on the MACHO photometry of 3674 RRO 
Lyrae stars toward the Galactic Bulge. Cook, Kunder, & Popowski (2007, in preparation). 
Cook et al. (2007, in preparation) divide the MACHO bulge RRO Lyrae variables into those 



1 RR0 stars have been tr aditionally cal l ed RR ab stars and are simply fundamental pulsators. For the new, 



intuitive nomenclature, see lAlcock et al.1 (|2000h 
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associated with the Galactic Bulge (3256 stars), and the ones associated with the neighboring 
Sagittarius dwarf galaxy (409 stars). Their classification is assumed in this work. 

The structure of this paper is the following. The method and resultant coefficient of 
selective extinction toward the Bulge is discussed and compared to other estimates in §2. 
Reddenings of the Bulge RRO Lyrae stars are determine d in §3 using their m inimum V- 
band light colors These reddenings are then compared to iPopowski et al.l (120031 ) reddening 
estimates and the outliers are discussed. 



2. Coefficient of Selective Extinction 

2.1. General approach 

As the bulge RR Lyrae stars span a much larger range in color than Sgr ones, a coefficient 
of selective extinction is determined using only this subsample. Assuming that the extinction 
affects the V magnitudes much more significantly than the distance spread, the intrinsic 
scatter in absolute magnitudes, and the scatter due to various observational errors, the 
slope of the stellar locus on a Vrr (mean V magnitudes of the RR Lyrae stars) versus 
(V — R)rr (mean color of the RR Lyrae stars) color-magnitude diagram (CMD) yields 



Rv,vr = AyjEiy — R) (IKunder et al.ll2006l ). In other words, if in the distance modulus, 



V RR -M v -A v = 5 log(d) (1) 

we assume My and d is approximately constant, then the apparent magnitude, Vrr, is 
proportional to Ay: 

Vrr ~ Ay ~ — — -. (2) 

twyR 

Taking the (V — R)q of each RRO Lyrae star to be approximately constant, 

Vrr ~ (3) 

*w,vr 

which is in the form of a general linear equation with Vrr as the independent variable, 
(V — R)rr as the dependent variable, and RyyR as the slope. It is shown in the following 
section that the above conditions are approximately met. 

Figure [D is a plot of mean Vrr, versus mean (V — R)rr. Artificial star tests from 



Alcock et al.l (119991 ) indicate that the limiting magnitude of the photometry in the V and R 
band is about 21.5, which is 0.5 magnitudes fainter than the faintest RR Lyrae star in this 
MACHO Bulge sample. The Monte Carlo analysis discussed later in this section indicates 
the limiting magnitude has little effect on this paper's analysis. 
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The sources contributing to the intrinsic dispersions of the average V and R apparent 



magni tudes are now discussed and quantified following the approach from IClementini et al. 



(120031 ) 



1. photometric errors. Since many MACHO fields partly overlap on the sky, some 
stars may be counted in the database twice. Using 184 double-represented stars, the internal 
photometric uncertainty in (V—R), V, and R was found by comparing colors and magnitudes 
evaluated in different fields. For the MACHO data toward the Galactic Bulge, the individual 
error in Vrr is ± 0.15 mag and the individual error in (V — R)rr is 0.04 (Cook et al. 2008, 
in preparation). § As the error in Vrr and (V — R)rr are correlated, upon multiplying the 
error in (V — R)rr, (0.04 mags), by the slope of the fit, (4.28), this value is added to the 
error in Vrr. The photometric error in the points in Figured] is thus 0.32 mags. 

2. absolute magnitude of RR Lyrae stars. It is generally assumed that the absolute 
magnitude of an RR Lyrae star has a linear dependence on [Fe/H], and that AMy(RR) 
I A[Fe/H] = 0.2 mag dex -1 . The dispersion around the mean value of [Fe/H] = -1.00 
dex found in Galactic Bulge RR Lyrae stars is oy = 0.29 mag. (Kunder et al. 2007, in 
preparation). This corresponds to an absolute magnitude dispersion of about 0.06 mag. 

The level of evolution off the zero-age horizontal branch (ZAHB) can also affect at RR 
Lyrae star's absolute magnitude. From the vertical height of the HB of a number of globular 
clusters of different metallicities, the dispersion around the average RR Lyrae luminosity 
due to the evolution off the ZAHB of e ach individual RR Lyrae is estimated as ~ 0.08 mag 



(IClementini et al.ll2003l : ISandagelll990l ) . 



The scatter in RR Lyrae absolute magnitude due to evolutionary phase, helium content, 
and alpha-element abundance [a/Fe], are less significant and are not accounted for. 

3. bulge distance spread. 

Using the standard Hernquist model density pro file with a scale length, cha racteristic 



density, and characteristic velocity of the bulge from IWidrow fc Dubinskil (120051 ) Table 2, 
the size of the bulge is ~ 1 kpc. As the galactic latitude and longitude of the MACHO bulge 
fields ranges from -1.5° to - 10° and 0° to 10°, respectively, a correction between the line of 
sight of the RR Lyrae stars in the bulge fields and the plane of the galaxy is negligible. A 
distance spread of 1 kpc corresponds to a AV^ of ~ 0.29 mags. 

Adding in quadrature all dispersion contributions, the scatter in V for a fixed (V — R) 



2 For more infor mation on the calibration of the MACHO Photometry RR Lyrae stars Database, see 
Alcoc k et al. especially pgs. 1547 and 1551. 



- 5 - 



is 0.44 mags, with the major source of this scatter resulting from the extent of the bulge 
along the line of sight. 

Fitting a straight line to the data must be done with caution, since errors are present in 
both independent and dependent variables and are furthermore correlated with each other. 
Ordinary least squares (OLS), which assumes no error in the independent variable, can 
therefore not be used here. Instead a straight line fit to the data is obtained using an OLS 
bisector fit to the individual data points. 

The OLS-bisector method has been used when both variables are subject to measure- 
ment error and it is not clear which variable should be treated as the independent and 
which as the dependent . It has been show n to outperform other approa ches in such cases 



(Ba bu &: Feigelsonlll992l ; llsobe et al.lll990l ). Using the SLOPES program (IFeigelson fc Babu 
19921 ). the data are fit using the OLS-bisector method. The slope found is 4.28 ± 0.04 and is 
shown as a solid line in Figured! Upon restricting the fit to those stars with (V—R) > 0.34, a 
slope of 4.34±0.04 is obtained. A further cut to remove the very red stars with (V—R) > 1.1, 
yields a slope of 4.49 ± 0.05. These fit values are given in Table [TJ the Rv,vr chosen as the 
optimal regression is Rv,vr = 4.28, as it accounts for all data values. Table [T] shows that the 
slope is quite sensitive to just a few points, although these formal errors for the OLS-bisector 
method are quite small. The error in Rv,vr of ±0.2 is chosen to encompass the range of 
slopes. 

The OLS-bisector fit to our data is further investigated with a simple Monte Carlo 
analysis. To duplicate the observed RR0 Lyrae data in this paper, 3500 E(V — R) values 
between and 1.5 are randomly selected. E(V — R) is converted to Av using a selective 
extinction coefficient, Rv,vr, which we vary. Assuming an RR Lyrae absolute magnitude 
of My = 0.6 and a distance to the Bulge of 8 kpc, an observed V magnitude is found. As 
the observed V magnitudes in the MACHO data have an error of ~ 0.5 mags, Gaussian 
random numbers are generated with a dispersion of 0.5 mags, and added to the simulated 
V magnitudes. Using (V — R) = E(V — R) + 0.28, and assuming the uncertainty in the 
measured (V — R) colors is 0.05 magnitudes, the simulated V magnitudes can be plotted as 
a function of simulated (V-R) magnitudes. We make sure this simulated plot looks similar 
to Figure [TJ 

Next a cut is made to remove all stars which have V and R magnitudes less than 21.5. 
The OLS-bisector method is then used to derive the slope of the fit before and after the 
magnitude cut. We find there is no significant change in the value of the slope after the cut. 
However, the slope found using the OLS-bisector method is not always the value of the slope 
that was input. This has implications on the derived Rv,vr value. Apparently the errors 
in the photometry in this paper are sufficient to throw off the OLS-bisector fit for a given 
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value. 

Figure [2] shows the difference in the true Rvyn value and the OLS-bisector slope, 
i?(v,y_R)-OLS slope, as a function of the true Rv,vr value. The scatter in this plot is mainly 
due to the errors in the ^-magnitudes and (V — R) colors. A slope with a value of ~ 4.6 
exhibits the least amount of deviation between the true and calculated Rv,vr value. The 
Rvyn value found in this paper is close to 4.6, so the slope found by the OLS-bisector method 
should not be significantly biased. Even values well above and below 4.6 do not differ from 
the true Rv,vr value by more than 0.2 mags. Other simulations with somewhat different 
input parameters (e.g., uncertainty in photometry) are performed and results similar, but 
not identical to Figure [2], are found. Thus, we can not simply correct the slope that is ob- 
tained from the real data using these simulations. Thus, an uncertainty of 0.2 in the Rvyn 
is adopted. 



2.2. Comparisons and Consistency 



It is an interesting question how this Rvyn value from the MACHO RR Lyrae stars 
compares with other reddening relations. Reddening relations determined from broad band 
photometric data depend on the intrinsic spectral energy distributions of the objects used 
to determine them as well as the particular filters. The MACHO and tm filters are 
considerably different from Johnson V and Kron-Cousins R. 



To compare our results to the standard extinction curve, the iFitzpatrickl (119991 ) ex- 
tinction curve with a constant total-to-select ive extinction parameter, Ry nv = 3.1, is used 
with the PHOENIX model synthetic spectra (IHauschildt. Allard fc Baronlll999l ). The John- 
son BVR transmission curves from Bessell (|l990|) yield R vyR = 4.88 for T eff = 6000 K, 
R VVR = 4.92 for T eff = 7000 K and R V y R = 5.07 for T eff = 30000 K. An effective temper- 
ature of 6000 K - 7000 K is typical of RR Lyrae stars, and an effective temperature of 30000 
K is typical of early type stars. The synthetic spectra adopted is one with [Fe/H] = — 1.0 
and logg = 2.5, but varying [Fe/H] and surface gravity does not significantly affect the Rvyn 
values. 

The MACHO bu and Tm filters are non-standard, but we re transformed into the Johnson 
V and Kron-Cousins R as described in lAlcock et al.l (Il999l ). We simulated this procedure 
by using the MACHO b M and r M transmission curves and our synthet ic spectra to obtain 
MACHO 6m and tm magnitudes. These were then converted using the lAlcock et al.l (Il999l ) 
transformation. It is found that Rvyn = 4.34 for T e ff = 6000 K, and Rvyn = 4.26 for 
T e ff = 7000 K. Both these values are in agreement with the value obtained in this analysis. 
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In comparison, a Rv,bv = 2.5 would yield Rv,vr = 3.77. 



Reddening laws are normally derived from differential measurement of reddened and 
unreddened luminous early-type stars. But such stars are absent in the Galactic Bulge. This 
analysis uses Population II type stars and indicates that the reddening law in the bulge is on 
average similar to the standard solar neighborhood value of Rv,bv = 3.1. Changing Rv,bv 
to e.g. a non-standard value, does significantly affect the measured Rv,vr- Quantitatively, 
it is found that the change in Rvvr is 68% of the change in Ry BV . 



Using Galactic bulge Red Clump Giants from the OGLE dataset. lUdalskil (120031 ) present 
substantial evidence that the ratio of total to selective absorption, Ryyi, is much smaller 
toward the Galactic Bulge than the value corresponding to the standard extinction curve 
and that Rv,bv varies considerably along different lines of sight. Their value of Ry,vi ~ 
2.1, as oppose d to the standard value of about 2.5, corresponds to an Rv,bv of about 2.6. 
Stanekl (119961 ) constructed extinction and reddening maps for Baade's window from color- 
magnitude diagrams obtained by OGLE which showed that extinction in Baade's window 
is quite i rregular, varying b etween 1.3 and 2.8 mag in Ay^vi with an estimated error of ~ 
0.1 mag. iRufHe et al.l (12004 ) use 70 Galactic planetary nebulae observed using narrow-band 
filters to find observed < Rv,bv > — 2.0 toward the bulge. 

The Rv,vr value determined from the RR Lyrae stars above, indicates that the standard 
Galactic reddening law can in general be adopted in studies of objects toward the bulge. We 
caution that the value of Rv,bv can vary in certain directions but that on average, the 
standard reddening law can be used. 



3.1. 



3. Reddening 
Determination of intrinsic (V — R)q colors 



Sturch 


(1966 


) and 


Blanco 


(1992) 



;ue that apparent (B — V) colors of RR Lyrae stars 
at phases close to the minimum light can be utilized to measure the amount of interstellar 
reddening along the line of sig ht to the star s ince t he intrinsic (B — V)q colors seem constant 
in the phase interval 0.5 - 0.8. iMateo et al.l (119951 ) suggest that (V — I) colors might also be 
used, perhaps with a smaller or negligible metallicity correction. We investigate this further 
for (V — R) colors below. 

The colors at minimum light, (V — R)q, can be computed via two different methods. In 
the first method, the Fourier fits to both V and i?-band data are used to average the colors, 
in magnitude units, for phases between 0.5 and 0.8. However, we find that when (V — R)q 
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is determined using this method, it correlated with amplitude. In the second m ethod, the 



(V — R) color is found at minimum ^/-band light. iKanbur fc Fernando! (12004 ) used this 
technique and found that intrinsic color at minimum ^/-band light is nearly independent 
of period and V^-band variation amplitude. Since this method yields (V — R) values that 
are very weakly correlated with amplitude, this is the method employed here. Table [2] 
presents our determination of unreddened (V — R)o colors for eleven well observed field RRO 
variables (the data in Table [2] come from Liu & Janes 1989; Cacciari et al. 1987; and Layden 
2005, private communication). The color excess, E(V — R), is computed by multiplying 
the published E(B — V) values by the ratio of Rvrv I Ryypt = 3.315/ 5.16, where the values 
R VBV = 3.315 and R vyR = 5.16 are taken from ISchlegel ' et all Jl998f )PI A plot of E(V - R) 
vs (V — R) has a linear character. This is quite remarkable, considering some of these stars 
exhibit the Blazhko effect (e.g., AW And), some are affected by both the Blazhko effect 
and a "phase-lag" problem (e.g., SS For), and some have evidence o f shock waves in their 



atmosp heres (e.g., RV Phe, V440 Sgr) as discussed by lCacciari et al.l (Il989a[ ). ICacciari et al. 



( 11989bl ). It has been shown that in many RR Lyrae stars the Blazhko effect does not affect 
colors at minimum light, making this method robust. We conclude that the mean (V — R)q 
color of RRO Lyrae stars at minimum V^-band light is 0.28 ± 0.02. Table [2] suggests that the 
intrinsic color at minimum light is independent or at least very insensitive to metallicity. 



3.2. Reddening toward the bulge RR Lyrae stars 

The observed and intrinsic colors at minimum V^-band are put forth to determine the 
color excess of the bulge RRO Lyrae stars: 

E(V-R) = (V-R)-(V-R) Q . (4) 

The Fourier fits performed to find (V—R) color at minimum V^-band light are sensitive to 
the quality of a lightcurve. To insure high E(V — R) accuracy without sacrificing the sample 
size, the Fourier fits in both passbands for all stars are individually examined. Only those 
with fits that approximated the RR Lyrae light curve nicely are taken (i.e., if the Fourier 
method mimicked features in the RR Lyrae lightcurve well, such as the dip and sharp rise 
of the curve at minimum light). This criterion results in 3525 RRO Lyrae stars with well 
determined reddening values. Example entries of this reddening catalog are presented in 
Table 3; the complete listing is available in the electronic version. 



3 We explicitly assume here that for this sample of RR Lyrae stars, extinction is on average well represented 
by the standard Rv,bv value. 



- 9- 



Multiplying the 3525 E(V — R)rr values by the selective extinction coefficient, Rv,vr, 
allows the determination of the visual extinction. A map of the visual extinction, Ay, is 



shown in Figure [3J As iPopowski et al.l (120031 ) also remark, it is immediately obvious that on 
large scales, extinction is regularly stratified parallel to the Galactic plane. 

To determine the distance modulus to the Galactic Bulge from the RRO Lyrae stars, 
an RR Lyrae absolute magn itude of M v = 0.59 ± 0.03 at [Fe/H] = —1.5 as compiled by 
Cacciari fc Clementinil (120031 ) is adopted. If a metallicity dependence of My = (0.25 ± 
0.05) [Fe/H] +constant is adopted (consistent with the range of slopes found in the literature), 
then using the average metallicity of [Fe/H] = —1 from above, an My = 0.72 ± 0.04 for RR 
Lyrae stars in the Bulge is obtained. This predicts (m — M)o = 14.7 ± 0.1 for Bulge RR 
Lyrae. Adopting t he RR Lyrae statistical p arallax calibration of My = 0.77 ± 0.13 at 
[Fe/H] = —1.6 from lGould fc Popowskil (119981 ). then using the same metallicity dependence, 
an (m — M)o = 14.5 ± 0.16 for RR Lyrae stars in the Galact ic Bulge is predic t ed. T his is 
identical to the Galactocentric distance modulus measured by lEisenhauer et al.l (120031 ) from 
the orbit of a star around the central black hole. 



3.3. Comparison with complementary MACHO determination 



In this subsection the reddenings obtained here are tested against IPopowski et al.l (120031 ) 
values, which are based on photometry from the MACHO bulge fields. The MACHO data of 
the Galact ic Bulge are based on o bservations of 94 bulge fields, each field covering an area of 
43' by 43'. IPopowski et al.l (120031 ) use 4' x 4' sky regions (tiles) as resolution elements. They 
show that mean colors for these tiles can be converted to extinction, and can thus be used 
to derive a visual extinction map. Each RRO Lyrae star is matched with the corresponding 
tile, an d the color exces s , EjV — R)jjR, of the line of sight to the RRO Lyrae sta r is compared 
to the Popowski et al. ( 2003 ). E(V — R)cmd, as calibrated to Stanek ( 1996 ) extinction^. 
Figure H] shows this linear relationship, with the line of slope unity over-plotted. The E(V — 
R)r r values become large r than E(V — R)cmd when reddening increases, which could mean 
the IPopowski et al.l (120031 ) extinction map underestim ates the extinction in areas with high 
reddening. Performing the iBessell fc Germany! (119991 ) calibration to the MACHO RR Lyrae 
lightcurve data does not alter this result. An increase in scatter is expected in stars with high 
reddenings, since generally a redder RRO Lyrae star is fainter and is more likely to have noisy 
photometry and lightcurves. Especially near the Galactic center, differential reddening due 



4 The Popowski et al. ( 2003 ). E(V — R)cmd, as cali brated to Dutra et al. ( 2003h extinction is about 0.07 
redder than E(V — R)cmd calibrated to lStanekl (|l996l ). 
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to nonuniform distribution of the intervening dust on small spatial scales makes the study of 
the s te llar population diffi cult even in the infrared (see, e.g.. iNarayanan et al.lll996l ; iDavidge 
19981 ) . iFrogel et al.l (119991 ) show in their Figure 5, that the amount of differential reddening 
in a field is directly proportional to the average r eddening: the more re ddening present in 
a field, the more patchy that reddening is. As the iPopowski et al.l (120031 ) extinction map is 
based on the average color of 4' x 4' sky regions, their average colors and hence reddenings 
would be susceptible to differential reddening, explaining the slight discrepancy between 
their reddenings and ours at higher reddenings. 

What could cause the systematic bias in t he E(V — -R)rr values to become larger than 
th e E(V — R)c,M r> as the reddening increases? IPopowski et ajj (120031 ) perform a calibration 
to IStanek! (119961 ) extinction and this calibration assumes a uniform Rv,bv value. This as- 



sumption does not hold, ho wever, since the extinction toward the Bulge has nonstandard 
properties (Popowski 2000). lUdalskil (120031 ) find that the difference between the measured 
line of sight Ry,vi, and the standard value of Rv,vi depends on color. Furthermore, they find 
that this difference is larger for redder V — I colors and practically negligible in the range 
where the OGLE data were calibrated by standards (V — I < 2). Similarly, the difference in 
the coefficient of selective extinction Rv,vr and the standard Ry yu would also increase with 
increasing V — R and hence with increasing E(V — R) cmd (since IPopowski et al.l (120031 ) uses 
the average color to determine E(V — R)cmd)- As the Popowski et al. ( 20031 ) E(V — R)cmd 
calibration neglects a change in Rv,vr with increasing color, at higher E(V — R)cmd the 
values will be underestimated. 



3.3.1. Outliers 

The most obvious outlier in Figure H] is the RR0 Lyrae star that is designated in the 
MACHO database as 175.30921.95. It has a negative E(V—R)rr, which is unphysical. It has 
a period, amplitude and light curve shape that confirm it is a fundamental mode pulsator. 
However, as is clear from its position on a color-magnitude diagram, wit h (V — R) = 0.21 



it is one of the bluest stars in the dataset (see Figure [p. ISturchl (119661 ) suggests that as 
many as 3% of the RR0 Lyrae stars may have an unusual blueness associated with them. 
The only stars with a color bluer than this star and with an R-band lightcurve suitable for 
reddening determination, are stars 125.23192.175 and 117.26079.4609, with (V — R) = 0.14 
and (V — R) = 0.19, respectively. The bluer of these two stars' E(V — R)rr deviates 
by 0.3 mag from E(V — R)cmd, while the other one by only 0.06 mag. Other very blue 
stars include 129.26623.1139 with (V — R) = 0.24 and 146.28412.88 with (V — R) — 0.25, 
but their E(V — R)rr is in agreement with E(V — R)cmd- One more characteristic that 
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distinguishes the blue outlying stars 175.30921.95 and 125.23192.175 from their blue non- 
outlying counterparts, is the fact that these stars have periods that lie in the range 0.469 - 
0.472 days, a range in which unstable light curves are common (IKinman fc Carrettalll992l ). 
Spectra of unusuall y blue s t ars in the period interval of . 469-0 .472 days, AR Her and BB 
Vir, were taken by ISturchl (119661 ) and iFernley fc Barnes! (119961 ) . It was found that these 
stars have systematically stronger hydrogen lines at minimum light than all the other stars, 
indicating a higher tempera t ure. The cause of the higher temperature, however, remains 
unknown. IFernley fc Barnes! (119961 ) suggested that BB Vir could be a RR0 Lyrae star that 
for some reason has evolved a long way beyond the fundamental blue edge without yet 
changing mode. The "hysteresis" effect is well-known in globular clusters and it is possible 
that stars 175.30921.95 and 125.23192.175 are also extreme examp les of it. Star s at the 
blue edge of the instability strip should have a larger amplitude (jSandage 1981 ). Stars 
175.30921.95 and 125.23192.175, the blue stars with deviating E(V - R) values, have about 
twice as large V amplitudes as the other two blue stars, 129.26623.1139 and 146.28412.88. It 
is also possible that the MACHO photometry of the star 175.30921.95 and 125.23192.175 is 
incorrect. Follow-up observations would help clarify this point. In conclusion, although rare, 
reddening values for RR0 Lyrae stars determined by finding the (V — R) color at minimum 
V-band light can be in error if the star is unusually blue. This can be corrected if the 
temperature of the RR0 Lyrae star is known. 



4. Conclusions 



A quantitative analysis of reddening toward the Galactic Bulge is performed. Assuming 
that the extinction affects the V magnitudes more significantly than the distance spread and 
the intrinsic scatter in absolute magnitudes, we interpret the slope of the stellar locus on the 
V versus (V — R) diagram as the selective extinction coefficient, Rv,vr = Ay/E(V — R). 
From the MACHO RR0 Lyrae sample toward the Galactic Bulge, we derive Rv,vr = 4.3±0.2. 
This corresponds to the average value observed in the solar neighborhood Rv,bv = 3.1 ±0.3, 
and indicates that the optical reddening law toward the bulge is consistent with the standard 
extinction law. We show that the mean minimum-light (V — R)q color of an RR0 Lyrae star 
is nearly constant at minimum V^-band light and equals 0.28 ±0.02, regardless of metallicity. 
Using this property, reddening values to 3525 RR0 Lyrae stars are found. Parameters derived 
for individual stars are presented in the electronic version. 
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Table 1: Slope (R v ,vr) of 3256 Bulge RRO stars plotted on a V versus (V — R) CMD. 



Fit 


Slope 


Error 


zero point 


Error 


OLS bisector, all stars 


4.28 


0.04 


14.37 


0.03 


OLS bisector, (V - R) > 0.34 


4.34 


0.05 


14.32 


0.03 


OLS bisector, (V - R) >0.34, (V - R) <1.11 


4.48 


0.05 


14.24 


0.03 



Table 2: Derivation of minimum-light (V — R) colors of field RRO Lyrae stars 



Star 


(V-R) 


E(V - R) 


(V - R) 


[Fe/H] 


V Amp 


SWAnd 


0.34 


0.035 


0.295 


-0.24 


0.89 


RRCet 


0.30 


0.019 


0.301 


-1.45 


0.78 


RXEri 


0.34 


0.053 


0.277 


-1.33 


0.87 


RRGem 


0.35 


0.054 


0.276 


-0.29 


1.15 


RRLeo 


0.31 


0.034 


0.266 


-1.60 


1.23 


TTLyn 


0.30 


0.014 


0.286 


-1.56 


0.70 


AVPeg 


0.34 


0.056 


0.274 


-0.08 


1.08 


TUUMa 


0.30 


0.018 


0.272 


-1.51 


1.02 


SSFor 


0.27 


0.012 


0.248 


-0.94 


1.23 


RVPhe 


0.30 


0.006 


0.284 


-1.69 


0.71 


V440Sgr 


0.33 


0.076 


0.254 


-1.40 


1.16 



Table 3: Fragment of the catalog of reddenings for 3525 RRO Lyrae stars in the MACHO 
Galactic bulge fields 



Star ID 


Ra 


Dec 


Gal I 


Gal b 


Period 


V mag 


R mag 


V amp 




E(V - R) 


101.20648.645 


18:04:31.80 


-27:33:55.8 


3.286 


-2.904 


0.590062 


17.58 


16.77 


0.57 


0.808 


0.53 


101.20649.500 


18:04:21.36 


-27:27:44.6 


3.357 


-2.820 


0.570140 


17.11 


16.46 


0.70 


0.729 


0.45 


101.20650.531 


18:04:25.32 


-27:26:16.1 


3.385 


-2.821 


0.557547 


17.29 


16.44 


0.81 


0.723 


0.44 


101.20650.965 


18:04:24.24 


-27:27:40.7 


3.363 


-2.829 


0.505202 


18.38 


17.57 


0.94 


0.829 


0.55 


101.20652.684 


18:04:23.16 


-27:17:18.2 


3.512 


-2.740 


0.570760 


17.07 


16.36 


0.25 


0.720 


0.44 


101.20653.526 


18:04:25.32 


-27:11:42.4 


3.597 


-2.702 


0.469922 


17.26 


16.53 


0.76 


0.688 


0.41 


101.20654.508 


18:04:30.72 


-27:11:14.6 


3.614 


-2.716 


0.552941 


17.13 


16.33 


0.88 


0.824 


0.54 


101.20654.938 


18:04:32.88 


-27:08:09.6 


3.663 


-2.697 


0.457945 


17.51 


16.78 


0.70 


0.778 


0.50 


101.20778.1172 


18:04:36.84 


-27:35:11.4 


3.276 


-2.930 


0.514020 


17.76 


16.95 


0.99 


0.817 


0.54 


101.20779.378 


18:04:41.16 


-27:28:17.4 


3.385 


-2.888 


0.633424 


17.06 


16.26 


0.53 


0.835 


0.55 
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Fig. 1. — Color-magnitude diagram for the 3256 MACHO Galactic Bulge RRO Lyrae stars. 
The best fit mean V versus mean (V — R) relation is shown by the solid line: Vrr = 
(4.28 ±0.0A)(V - R)rr + (14.37 ± 0.03). and was obtained using the OLS-bisector method 
jFeigelson & Babul Il992h .. 
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Fig. 2. — The difference in the true Rv,vr value and the OLS-bisector slope, R(v,vr) - 
OLS-fitted slope, as a function of the true Rv,vr value. 
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Fig. 3. — Reddening map of A v values of the central Galactic region. The scale is given at 
the bottom. Insert presents three disk fields separated from the others by 10°. 
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Fig. 4. — A comparison of E(V — R) derived using (V — R) at minimum lig ht, E(V — -R)rr, 



for 26 05 RR Lyrae stars with color-magnitude based E(V — R)cmd from iPopowski et al. 



( 120031 ). The relation is linear, with E(V — R)rr values slightly larger than E(V — R)cmd at 
high extinction. 



